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I. INTRODUCTION 



Understanding the phase-space transport of magnetically-trapped relativistic electrons 
is an intrinsically interesting general problem, and it is an especially important problem 
in magnetospheric plasma physics because these energetic particles can damage spacecraft 
electronics and they present a radiation hazard to astronauts. The quasilinear theory of 
such anomalous transport processes is based on resonant wave-particle interactions in which 
characteristic wave frequencies match one or more of the orbital frequencies associated 
with the gyration, bounce, and drift motion of charged particles [TJ EJ E|- In particular, 
diffusion due to drift-resonances with low-frequency magnetohydrodynamic (MHD) waves 
[U El El and diffusion due to cyclotron-resonances with various high-frequency waves [7J El 
El CD] are frequently cited as important transport mechanisms for relativistic electrons in 
Earth's magnetosphere. Previous work on these transport mechanisms has typically been 
restricted to either the low-frequency wave interactions, which may break the second and 
third invariants, or the high-frequency wave interactions, which may break the first and 
second invariants. However, there has been little work on a unified theory in which all the 
interactions can be considered together. 

Our previous work (henceforth refered to as Paper I) investigated relativistic quasi- 
linear diffusion transport due to low-frequency electromagnetic fluctuations, which pre- 
served the first adiabatic invariant, based on the low- frequency relativistic gyrokinetic 
Vlasov equation [T2~l IT51 ITlj . In the present work, we modify the low-frequency gyroki- 
netic formalism to allow for arbitrary-frequency electromagnetic fluctuations (and, thus, 
allow for processes which break any one of the three invariants) and retain full finite- 
Larmor-radius (FLR) effects associated with the electromagnetic fluctuations. Although 
FLR effects are not important for interactions of relativistic electrons with MHD waves, 
they are crucial for interactions with cyclotron-frequency waves E] • 

A. Canonical Quasilinear Formalism 

The derivation of a quasilinear diffusion equation in canonical action space (with coordi- 
nates J) was first performed by Kaufman in which the background Vlasov distribution 
Fq(3,t) evolves on a slow time scale (r = e 2 t) determined by the amplitude (ordered by 
the dimensionless parameter e) of the electric and magnetic fluctuations. The quasilinear 
diffusion equation in action space is given in its general form as 



where the action coordinates J may either be exact invariants (corresponding to exact sym- 
metries) or adiabatic invariants (corresponding to approximate symmetries). The canonical 
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quasilinear (CQL) diffusion tensor [TJ |TTj is expressed as 



Dcql = Yj mm S(uj k -m-n) \5H mk \ 2 \ , (2) 

m,A; 

where 5H m k denotes the Fourier component of the perturbation Hamiltonian 5H with a 
discrete frequency spectrum (represented by the wave frequency u) k ) and m ■ fl = m' fli 
(summation over repeated indices is, henceforth, implied), with f2j = dH /dJ l denoting 
the unperturbed orbital frequency derived from a Hamiltonian Hq averaged over all fast 
orbital time scales. 

The expression (j2J) clearly establishes the paradigm of quasilinear transport theory: 
wave-particle resonances (where u k = m ■ f2) introduce explicit violations of the invariance 
(exact or adiabatic) of action coordinates leading to stochastic transport in action space. 
Note that the canonical quasilinear diffusion tensor (J2J) can be written 

n _ V h( cql ) f( cql ) (i\ 

L>CQL - 2^ d m,fc L m,k > \ 6 ) 
m,k 

in terms of the canonical quasilinear coefficient matrix = mm and the canonical 

quasilinear potential = it 5(ujk — m • ft) \5H m k\ 2 - Here, the dimensionality of quasi- 

linear transport is represented by the dimensionality of the coefficient matrix d mj fc and 
the universality of quasilinear transport is represented by the fact that a single quasilinear 
potential T m k describes quasilinear transport along different action coordinates. 



B. Organization of Paper 



The remainder of the paper is organized as follows. In Sec. II, we introduce the unper- 
turbed and perturbed relativistic guiding-center equations in terms of magnetic coordinates 
(ip, if, s) associated with an unperturbed axisymmetric magnetic field configuration. The 
invariant coordinates I = (J g ,s, Jd) in axisymmetric magnetic geometry are the relativis- 
tic guiding-center gyro-action J g , the relativistic guiding-center kinetic energy e, and the 
guiding-center drift action J d = qip/c. Here, energy e is used instead of the bounce action 
Jb to ensure that all three of the invariant coordinates are local (i.e., they depend on local 
properties of the background plasma). The perturbation guiding-center Hamiltonian 5H, 
expressed in terms of the perturbation electromagnetic potentials (5(f), SA), causes the de- 
struction of the unperturbed invariants I (i.e., 51 = {I, 5H} ^ 0, where { , } represents 
the Poisson bracket) and leads to stochastic quasilinear diffusion in invariant I-space due 
to resonant wave-particle resonances. 

In Sees. Ill and IV, an explicit expression for the relativistic quasilinear diffusion equa- 
tion Q in axisymmetric magnetic geometry in terms of the invariant coordinates ( J g , e, Jd) 
is derived. In Sec. Ill, the quasilinear coefficient matrix d m k is derived and then, in Sec. IV, 
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the quasilinear potential r m ^. The present derivation follows closely the derivation found 
in Paper I [TTJ. We summarize our results in Sec. V and discuss applications. 

Since quasilinear diffusion is often described in the literature (e.g., see Ref. [2]) in terms 
of the equatorial pitch angle instead of the gyro-action, the corresponding relativistic quasi- 
linear diffusion equation is presented in Appendix A. In particular, we show how the 
quasilinear coefficient matrix transforms when the relativistic guiding-center gyro-action J g 
is replaced with the equatorial pitch angle. Lastly, Appendix B presents two expressions for 
the bounce-averaged drift frequency which might be useful for computational applications. 



II. RELATIVISTIC GUIDING-CENTER DYNAMICS IN AXI SYMMETRIC MAGNETIC 
GEOMETRY 



We begin our analysis of relativistic quasilinear transport in axisymmetric magnetic 
geometry by first providing a general representation of axisymmetric magnetic geometry 
in terms of magnetic coordinates (also known as Euler or Clebsch potentials). Next, we 
present the equations of relativistic guiding-center Hamiltonian dynamics in unperturbed 
and perturbed axisymmetric magnetic geometry. Unperturbed guiding-center Hamilto- 
nian dynamics establishes the existence of three invariants, which are then destroyed by 
arbitrary-frequency electromagnetic perturbation fields. 

In the spirit of perturbation theory, we assume that electric fields are not part of the 
quasi-static background fields (i.e., electric fields are automatically viewed as perturbation 
fields) and any axisymmetry-breaking magnetic fields are viewed as zero-frequency magnetic 
perturbation fields. The validity of these assumptions will be checked in future work when 
an alternate model including a more realistic non-axisymmetric magnetic geometry and a 
background electric field is considered. 

A. Axisymmetric Magnetic Geometry 

The background magnetic field in axisymmetric magnetic geometry can be expressed 
in terms of magnetic coordinates (ip, <p, s) as 



where ip denotes the magnetic (radial) flux function (i.e., B ■ Vt/> = 0), the azimuthal angle 
<p is an ignorable coordinate (i.e., in axisymmetric magnetic geometry, unperturbed scalar 
fields are independent of (p) , and s is the parallel spatial coordinate along a single magnetic 
field line labeled by (ip,(p). 

The magnitude B(ip,s) of the magnetic field, on the one hand, is defined from Eq. (0J) 
as B = V"0 X Vy? ■ Vs and i? _1 (^, s) denotes the Jacobian for the transformation X — > 




(4) 
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(if),(p,s). The local magnetic unit vector b = B/Z? = <9 S X, on the other hand, can be 
expressed as 

b = Vs + a(ip,s) Vip (5) 

Here the scalar field 



a(ip,s) = 




|V^| 2 



which characterizes the non-orthogonality of the coordinates (ip,s), is required to ensure 
that Vxb^O and is associated with magnetic curvature [TT] through the relation d s b = 
(d s a) Vip. We note that, for an axisymmetric dipole magnetic field (which is curl-free), we 
can write B = Vx(/0, s) and, thus, using Eq. ©, we find dx/ds = B and dx/dip = aB. 

Next, we define the radius R(tp,s) = |Vv?| _1 (which measures the distance to the 
symmetry axis) and obtain the expression |V^| = BR from Eq. (J1J). We can, thus, define 
the following perpendicular unit vectors 

$ = (BR)- 1 Vip and (p = R Vp, (6) 

such that the unit vectors (ip, (p, b) form a right-handed set of unit vectors (i.e., b = if) X (p). 

B. Unperturbed Guiding-Center Dynamics 

The unperturbed relativistic guiding-center dynamics of a charged particle of mass M 
and charge q in axisymmetric magnetic geometry is represented in terms of the magnetic 
coordinates (ip,<p,s), the unperturbed relativistic guiding-center kinetic energy 

e = ( 7 -l)Mc 2 = ^M 2 c 4 + 2 Mc 2 J g u g + p\c 2 - Mc 2 , (7) 

where u g = qB/Mc denotes the rest-mass gyro-frequency and p\\ = j Mv\\ denotes the par- 
allel component of the relativistic kinetic momentum, the relativistic guiding-center gyro- 
action J g = |p_i_| 2 /(2 Muj g ), and the gyro-angle (. The unperturbed relativistic guiding- 
center Lagrangian is written as 

L = J d p + p\\ s + J g C - e, (8) 

where A = ip Vy? is the vector potential associated with the magnetic field (@J) and Jd = 
qip/c denotes the drift action, while, using Eq. (J7J), the magnitude of the relativistic guiding- 
center parallel momentum is 

hi = 7 M|«||| = ^e 2 /c 2 + 2M(e - J g u g ). (9) 

The unperturbed relativistic guiding-center Poisson bracket { , }, on the other hand, is 
expressed in terms of two arbitrary functions F and G on guiding-center phase space as 

rrn ^ dFfdG dG\ fdF dF\ dG 
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dF_ (dG_ dG ca dG\ ( 3F_ dF_ ca dF\ dG 

dpi \dJ d U ' d de q ds j \dJ d Ud de q ds J dp 

fdF dG dF_ dG\ 

" \ds de de ds J ' 

where uj c = 7 _1 u g denotes the relativistic gyro-frequency and the azimuthal drift frequency 
is defined as 

/ du c ac du c \ V\\ da 

The Jacobian for the guiding-center transformation 

(x,p) -»• Z a = (J d ,tp, s,e, J g ,C) 

is J = l/|f|||, where \v\\\ is obtained from Eq. ©. We note that the Poisson bracket (fTU|) 
can also be expressed in phase-space divergence form as 

jM jfj *&)- (12) 

where J af3 = {Z a , Z 13 } denotes the elements of the antisymmetric Poisson-bracket tensor 
and we used the Liouville identities (f3 = 1, 6) 

These identities follow from the divergenceless property of the Hamiltonian flow (associated 
with Liouville's Theorem): 

d ■ 1 d , „ . x 1 d ( „ mtt dH\ 



dZ J dZ a V > J dZ a \ dZP 



J dZ a ^ > dZ? dZ a dZ?' 
Since the Poisson matrix J a @ is anti- symmetric, the term J a @ d^pH vanishes, and we find 

o = ^^(JJ^ dH 



J dZ a v / dZP' 

which must be true for all Hamiltonians H and, thus, we obtain the Liouville identities 
(I13|) . In particular, the Liouville identities for p> and (, respectively, are 

d f 1 \ d ( ca \ + d_{ q 



dJ d \ \v\\\J ds \q\v\\\J de 



(14) 



d ( 1 \ d ( uj c 



dJe \ \v\\ \ I de 



which follow from the definitions for u c , \v\\\, and We further note that although the 
magnetic coordinate ip and the drift action Jd are simply related in axisymmetric magnetic 
geometry, this simple relation is destroyed by non-axisymmetry; we, henceforth, use ip 
and Jd separately as follows: Jd is used whenever dynamics is concerned while ip is used 
whenever magnetic coordinates are concerned. 

As in Paper I [11], the unperturbed guiding-center Hamilton's equations, expressed 
in terms of the coordinates ( Jd, (p, s; e, J g , () and the unperturbed relativistic guiding- 
center Hamiltonian Hq = e, include the magnetic-coordinate components of the relativistic 
guiding-center velocity 

j d0 = 0, Lpo = u d , and s = v\\ (15) 

and the velocity-space guiding-center equations 

e = = j g o and Co = • (16) 

Thus, the drift action Jj (= qip/c), the kinetic energy e, and the gyro-action J g are the 
three invariants of unperturbed relativistic guiding-center motion in axisymmetric magnetic 
geometry. From these unperturbed Hamilton's equations we construct the unperturbed 
time evolution (Vlasov) operator 

d _ d d d d 

Jt = m +v "dS + UJd d^ + UJc d-c (17) 

In the absence of electromagnetic fluctuations, the unperturbed Vlasov equation d F /dt = 
0, therefore, implies that the unperturbed (background) Vlasov distribution F (I) is a 
function of the three guiding-center invariants I = (J g , e, Jd) associated with axisymmetric 
magnetic geometry. In the present work, in order to correctly account for wave-particle 
gyroresonance effects, we retain full finite-Larmor-radius (FLR) effects associated with 
electromagnetic fluctuations and, thus, we need an expression for the gyroradius p. The 
gyroradius vector p(X, s, J g ,C) is defined locally in terms of the magnetic coordinates as 



p = p (cos C ip — sin C, (p 



where the unit vectors (ip, (p) are defined in Eq. © and the magnitude p = v±/uj c of the 
gyro-radius vector p is given as 



= imjhy (19) 

In Sec. IV, the spatial dependence of the perturbed electromagnetic fields (in guiding- 
center phase- space) will involve the combination ip 1 + p- V"0 l , which denotes the position 
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x = X + p of a particle in terms of the guiding-center magnetic coordinates ip l = (ip, ip, s). 
Here, using Eqs. (JBJ) and ()18|). we find 

p-Vip = p|V^| cosC = A^/j cosC 



p-V<p = -p|V<p| sinC = ~(p/R) sinC 

p • Vs = — a p • = — a A-0 cos £, 
so that a generic perturbation scalar field 

<5x(X + p,t) = exp(p- V)5x(X,t) = 6gcx(il) } <p } s) e,J E ,(,t) 
can be expressed as 



(20) 



<^x(X + p,t) = exp 



ds 



(21) 

6 X (^,(p,8,t). (22) 



This expression will be used later in obtaining Fourier-Bessel expansions in gyrorangle £ 
for the electromagnetic perturbation fields. 

C. Perturbed Guiding-Center Dynamics 



The presence of fluctuating electric and magnetic fields 

SE = - V5<f) - cT 1 d t SA and SB = V x 5 A, 

where S<p and 5A are the perturbed scalar and vector potentials, respectively, implies that 
the three unperturbed invariants ( J g , e, Jd) are no longer invariants. The perturbed electro- 
magnetic potentials introduce perturbations in the relativistic guiding-center Lagrangian 
©: Lq — > L + SL, which, to first order in the perturbed potentials, yields 

Q 



SL 



S gc A- (X + p 



q5, 



(23) 



where the notation S gc x is defined in Eq. (|21J) . As a result of the magnetic perturbation (the 
first term in Eq. (J2BJ)), the relativistic guiding- center Poisson bracket (fTUJ) is also perturbed 
unless we define the perturbed Hamiltonian SH as 



SH = -SL 



q <L C 



S BP A \\ 



q 



(24) 



The destruction of the three unperturbed guiding-center invariants is expressed in terms of 
the perturbed guiding-center Hamiltonian SH and the unperturbed Poisson bracket (fTUj) as 

dSH 



SJ S 


= R, SH} 


Si 


= {e, SH} 


SJ d 


= {Jd, SH} 



I " Jt) 5H > 
__ dSH 

dip 



(25) 
(26) 
(27) 
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where the unperturbed time evolution operator d /dt is defined by Eq. (jT7jl . while the 
remaining perturbed guiding-center Hamilton's equations are 




<SC = {C, SH} = — + o>„ ^- «f, (28) 



ft = {., «*} = (v, » + = £) (29) 

* - OT > = " 7 S + <" I) w ' (30) 

From Eqs. (J23J) and (|2T|). we note that the gyro-action J g and the drift-action J d are de- 
stroyed if the perturbation Hamiltonian SH depends on the gyro-angle ( and the azimuthal 
angle (p, respectively. Moreover, the destruction of the gyro-action invariance leads to the 
possible loss of trapped particles while the destruction of the drift-action invariance leads 
to radial transport. 



III. RELATIVISTIC QUASILINEAR DIFFUSION TENSOR 



In this Section, we proceed with a two-time-scale analysis of the perturbed Vlasov 
equation 

e> d ^ + e d -f^ = -e{(F + e5F),5H h (31) 

where the unperturbed evolution operator do/dt is defined in Eq. (IT7j) and the Vlasov dis- 
tribution F is decomposed as F = F (I,t = e 2 t) + e5F. The fast-time-scale evolution 
equation shows, on the one hand, how the perturbed Vlasov distribution 5F evolves under 
the influence of electromagnetic perturbation fields. The slow-time-scale evolution equa- 
tion, on the other hand, shows how the background distribution Fo changes as a result of 
wave-particle resonances. Following our previous work (Paper I) , we now proceed with the 
separation of the perturbed guiding-center Vlasov distribution 

BF 

5F = —^5H + 5G (32) 
oe 

in terms of the adiabatic part of 5F (the first term) and the non-adiabatic part of 5F 
(the second term), which explicitly represents resonant wave-particle effects. The fast-time 
evolution equation for the non-adiabatic part 5G is obtained from the linearized Vlasov 
equation d 5F/dt = — 51 • dF /dl as 

d 5G _ d 5F dF d 5H 



dt dt de dt 

d dF d dF d 




dJd dip dJ g d( de dt 



8H = f SH, (33) 
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where the perturbed Hamilton's equations for SI = (SJ g , Si, SJ^) are given in Eqs. ([25 )1 - (|27 |) . 
and we have defined the operator T . 



A. Slow-Time Evolution 



Since the guiding-center background distribution F is quasi-static on the wave time 
scale and is independent of the azimuthal angle p and the gyro-angle (, we introduce an 
averaging operation [denoted as (• • •) and refered to as wave averaging] with respect to the 
fast wave-time-scale, the azimuthal angle, and the gyro-angle, with the property 



F = F„(I; r = e 2 t). 
Applying the wave-averaging procedure on Eq. (|31|). we obtain 



dFo 
dr 



({SF, SH}) 




j(SF SI 



d_ 

ds 



J(SF 6s) 



(34) 



(35) 



where the perturbed Vlasov distribution SF and the perturbation Hamiltonian SH have 
zero wave-averages: SF = = SH and expressions for SF = (SJ g , Si, SJd) and 5s are 
found in Eqs. (j25|l - (|27)l and pi. respectively 

When the slow-time evolution equation (|3"K|) for the guiding-center background distribu- 
tion Fq is expressed in terms of the non-adiabatic part SG of the perturbed guiding-center 
distribution, we find (after some algebra) 



dFh 
dr 



d 



SG 



dSH \ 



d 



dSH\ 




+ 



d_ 




+ 



d_ 

ds 



ac 



q\v\\ 



SG 



dSH \ 




^ dSH\ ( s „ dSH\ 



dSH\ o d 2 F 



de ) 



2 de 2 



(SH 2 



(36) 



where a - 
identities 



= ± 1 and we made use of the Liouville identities (jl4)l as well as the 



(SFd^SH) = (SGd ( SH) and (SFd^SH) = (SGd^SH). 

The slow-time evolution equation (j3T?j) for the background distribution F contains terms 
associated with exact derivatives in invariant I-space (as expected) plus a term involving 
parallel spatial gradients (b- V = d/ds) along magnetic field lines. Since the background 
distribution F is independent of the parallel spatial coordinate s (i.e., dF /ds = 0), we 
must remove the parallel-gradient terms on the right side of Eq. (J36|) by introducing a 
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second averaging operation. 



B. Bounce Averaging and Fourier Decomposition 



To remove the parallel-gradient terms on the right side of Eq. (|36p. we introduce the 
bounce averaging operation 



ds J 



(37) 



where J7" = l/|f||| is the Jacobian introduced above, sz(I) and st/(I) are the turning points 
where the trapped particle's parallel velocity v \\ vanishes, J2a denotes a sum over the two 
possible signs of vn = ± \v\\\, and Tb denotes the bounce period 



7b 



E 



J ds. 



The bounce-average operation defined here yields the following identities 



--[J 
,J ds [J 



I 1 d 

and (jW^ J 



1 A 



7b 



Hence, by bounce averaging the slow-time evolution equation (jHEJ) while using these iden- 
tities, we obtain 

dF 1 d 



Or 



- W U(SGSI< 



(38) 



where we note that (5P) = and the bounce period Tb now appears as the new Jacobian. 

We now introduce the Fourier decomposition (assuming a discrete frequency spectrum 
{u>k} for the waves) 



/ 5G 
V 5H 



E E 

k m,i=—oo 



SG mek (s,a; I) 
5H m£k (s,a; I) 



exp z (rrup + £( — u) k t) 



(39) 



so that the fast-time evolution equation (J33|) for the non-adiabatic part SG can be written 

as 

d 



V " ds ~ \ UJk ~ m Ud ~ c 



oG m ch = £ 5G 



inlk 



iT 8K 



miki 



(40) 



where the differential operator T [defined in Eq. (|33jl ] becomes 



. OF, dF dF 

% I 771 "tttt - + « TT-r + ^fc 



9J d 



<9J, 



g 



ds 



(41) 
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Substituting the Fourier decomposition into Eq. (|38p. with Eqs. (J25J) and (|27)l. respec- 
tively, we find 



5G5J? 



- £Im($G mek 8H* mlk/ . 

m,l,k 



m Im {5G mlk bH* mik 

m,£,k 



5G 5Jd 

After performing various integrations by parts using Eqs. (j2*Uj) and ()33|). we also find 



(42) 
(43) 



SG Si 



d n SG dSG 



dt 



dt 



5H 



F5H - 



dSG 



5H 



05G 



5H 



u k Im iSGmik 5H* 



mlk / ' 



(44) 



where we used the identity ((F5H) 5H) = 0. We note that Eqs. P^l - fP^ all involve sums 
containing the term lm(5G m £k 511^^), so following our earlier work (Paper I) we introduce 
the quasilinear potential 



J- mlk 

so that we may replace 

Im (bGrnlk ^^mik 

into Eqs. @2J-©- 



T- 1 Im (5G mIk SHI 



mlk urL m£k/ i 



f dF dF «9iV 



(45) 



(46) 



C. Relativistic Quasilinear Diffusion Tensor 



By substituting Eqs. 
linear diffusion equation 



and (J4l)j) into Eq. (J3~8~j) . we obtain the relativistic quasi- 



dF (I,r) 
Or 



r h dF 



dF (I, T y 
dF 



(47) 



where the invariant coordinates are the gyro-action F = J g , the particle guiding-center 
kinetic energy I 2 = e, and the drift action J 3 = In Eq. (JUJ), the relativistic quasilinear 
diffusion tensor Dql has the following symmetric form 



D 



QL 



E 

m,l,k 



f i 2 £tu k im 
u k i uj\ ijj k m 
\ mi muo k m? 



mlk' 



(48) 
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Note the simplicity of the quasilinear coefficient matrix d mk when expressed in terms of the 
invariant coordinates I = (J g ,e, Jd). Since the quasilinear transport of trapped particles 
(e.g., see Ref. 0) is often discussed in terms of the equatorial pitch angle instead of the 
gyro- action J g , an alternative representation of the quasilinear coefficient matrix expressed 
in terms of the equatorial pitch-angle instead of the gyro-action is presented in Appendix A. 
Although the pitch-angle formulation will facilitate comparison with earlier works, however, 
the simplicity of the quasilinear coefficient matrix in Eq. (J48j) is lost. 



IV. QUASILINEAR POTENTIAL 



To complete the derivation of the relativistic quasilinear diffusion tensor Dql, with com- 
ponents given by Eq. (|48|). we must now solve for an explicit expression for the quasilinear 
potential T mm . 



A. Guiding-Center Perturbation Hamiltonian 



In the present work (as in Paper I [Uj), we use the gauge condition 8 A • Vy? = so 
that the perturbed vector potential is written as 

5 A = <L4|| Vs - 8(3 W, (49) 

where 8 An = d5a/ds denotes the parallel component of the perturbed vector potential and 
the perpendicular components of the perturbed vector potential 8A± = 5 A — 8Au b are 

8A ± = - \ 8fi + a Saj W. (50) 
Hence, the parallel component of the perturbed magnetic field is 

while the parallel component of the perturbed electric field is 

r L rj 1 88 A\ 35$ 

where 5$ = difi + c^dtba denotes the effective perturbation scalar potential. Note that the 
Fourier components of the perturbed parallel electric and magnetic fields are given as 



8E\\mk = 7^— = - [8(pmk ~ i — 8a mk ) , (53) 



8B\\ mk = imB (sPmk + a ~^~]- ( 54 ) 
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As in Paper I [TT], it is possible to find expressions for the perturbation potentials S(f> m k, 
Sa m k, and 5f3 m k in terms of electric covariant components 5Ei m k = <5E mfc • <9jX and magnetic 
contravariant components 5B l mk = <5B m fc • V?/>\ 



B. Fourier-Bessel Expansions 

Using the gauge condition introduced in the previous Section, the perturbation Hamil- 
tonian (J24j) becomes 

5H = q ^5 gc (j) - 5 gc aj + ^ v ± • (^ gc (3 + a 5 gc o^j . (55) 

We now consider the perturbation potential 

8 X (X + p,t) = J2 e'^SxkiX + p), 

k 

which, using Eqs. (f?T]) - (f22|) . yields the azimuthal- angle Fourier expansion 

oo 

5 Xfc (X + p) = e im ^[exp(-^sinC + cosCA) s)], (56) 

m=— oo 

where 

v ee rat and A = _ a . (57) 

We note that the argument A is a differential operator acting only on the perturbation 
fields. 

Next, we introduce the following Bessel-function identities 

oo 

exp (-i r) sinC) = E exp(- ijQ J 3 (rj) (58) 

j=-oo 
oo 

exp(cosCA) = J2 ex P(t/'C) IfW, (59) 
y=-oo 

where J 3 and Jy denote the Bessel and modified Bessel functions of order j and /, respec- 
tively, and 

i„ I oo \ / oo \ oo 

-/ e"^ £ e-*J,fo) £ ^/y(A) rfC = E 4^- (60) 

Z7F J V=-oo / \/=-oo / J=-oo 

The Fourier gyroangle expansion of the generic perturbed potential ()56j) is, therefore, ex- 
pressed as 

oo I oo \ 

5 Xfc (X + p) = £ e^+^ 2 J^) / J+/ (A) <W0M 

m,£=— oo yj=— oo y 

oo 

= E e^"^^, (61) 



m, l=—oo 
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where the Fourier-Bessel components Sx m ek( s ] I) are functions of the parallel spatial coor- 
dinate s and the invariant coordinates I = (J g ,£, Jd)- 

Lastly, using the expression (|54~jl for 5B\\ m £ k , the Fourier-Bessel expansion of the third 
term vj_ • 5A± in the perturbation Hamiltonian ([55)1. with 5A± given by Eq. ([50)1 . is pro- 
gressively transformed as follows. First, we begin with 



V ± -5A_L fc 

c 



( 9 

| V'0| v± sin ( \S/3 k + a — Sa 



ds 



■Rp± sinC E 



m 



(E J 3 tj+e 



6B\ 



mk 



B ' 



where Eqs. (I54|) and (|6*T|) were used, with the definitions IV^I = BR and p± = j Mv_\ 
Next, by substituting sin( = (e^ — e~^)/2i, we obtain 

q pi(m<P+ZC) 

c 



q ^ 

- v ± • 5A ±k = u c Rp ± 



ml 



2m 



5B 



\\mk 



B 



where the j-summation was re-arranged. Lastly, we use the Bessel recurrence relation 
Jj-iiv) ~ Jj+iiv) = 2rf x Jj(t]) to obtain 

" 6B n 



v_l • 5A ±k 

c 



m,£ 



E i*t$) ( - - J'M 



V 



>\\mk 



B 



oo r jd 



m, t=— oo 



In summary, the Fourier-Bessel components of the perturbation Hamiltonian (|55|) are 

5B\\ m £k 



hi 




c 


ds J 



+ Jjjj, 



gU/ C 



B 



where 



84> m ek \ 


oo 




= E 


Sa m ek J 


j=-oo 




oo 


5B\\ mtk 


= E 




j=-oo 



I S<p mk (ip,s) 
\ 5a mk (tp,s) 



-j;(v)) I J+£ (X)6B llmk (iP,s) 



(62) 

(63) 

(64) 
(65) 



The Fourier-Bessel components 5H m £ k (s,a; I) of the perturbation guiding-center Hamilto- 
nian can now be used explicitly in the expression for the quasilinear potential (|45|). 



C. Solution of Fast-Time Evolution Equation 



Following an approach detailed in Paper I, we remove the cx-dependence of 5H m g k (s, a; I) 
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by introducing a new perturbation Hamiltonian 

&K mtk = 5H mfk + - £ 5a m ek 

c 

= q5^ m ek + i-^ (jnuj A + iu^ 5a m£k + J s u c SB ^ ek (66) 
where the differential operator £ is defined in Eq. and 

S^ m ek = Hmik - i (ujk/c) 5a m ek = - 5Eu mek ds. 



Note that, as in Paper I (Appendix B), the new perturbation Hamiltonian can also be 
expressed as 



5K mek = — (muj d + £uj c - uj k ) [ SB ™ ek ds + J g u c SB \[™ £k + ilsEpmtk, (67) 
mc J B B m 

in terms of components of the perturbed electric and magnetic fields. 

We also introduce the new nonadiabatic part SG' m£k defined as 

SG' m£k = 5G m ik + i-Tba mik (68) 

c 

so that the fast time scale evolution for SG' m£k is obtained from Eq. as 
£ &G' mlk = £ 5G mek + i - T £ 5a m i k 



/ ~ q 



= iTySHmjik + i-£5a mek j = i J 1 ' 5K mek , (69) 

where we used the fact that the differential operator £ commutes with T . 

The rest of the analysis leading to the solution of SG' m£k in terms of the modified per- 
turbation Hamiltonian SK m £ k follows Paper I [TT]. First, we note that the transformation 
(5G m e k ,5H m £ k ) — > (5G' mek ,5K m £ k ) leaves the quasilinear potential invariant: 

T mek = T~ x Im (5G' m£k 5K* m£k ) , (70) 

as follows from properties of the averaging operations ((•■■)) an d the differential operator 
£. From Eq. (61) of Paper I, where we now replace moj^ — > {moj A + £u c ), the solution of 
Eq. is written as 



! J St, |f||| 

cot Q^SKmtk cosO^ + (dKmik sinfl)) \, (71) 
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7b 



where the s-dependent angle 9(s) is defined as 



ds' 



9(s; I) = / 1 — - u k - mu d (s) - £uj c {s') 



while the s-independent angle is defined as 



6(1) = u k - £ (uj c ) - m (u d ) ) . 



(72) 



(73) 



Here, (u c ) and (uj d ) denote the bounce-averages of the gyro-frequency and drift-frequency, 
respectively; in Appendix B, we present a simple expression for the bounce-averaged drift 
frequency (u d ) in terms of the bounce-averaged radial gradient (du c /dJ d ). 

Using simple relations derived in Paper I, we ultimately find the quasilinear potential 
^mik given as 



r m <? fc (I) = y (6K m ek{s; I) cos9(s; 1)^ Im(-cot0). 



Here, we note that 



(74) 



y cot 6 = ]T (u k -£ (u c ) - m (u d ) -nu h 



where u\ } = 2-n/Ty, denotes the bounce frequency. 
Lastly, using the Plemelj formula, we find 

— Im(-cote) = n5(L} k - m-(a> m }) 

n=— oo 

where m = (£, n, m) and 

m ■ (u> m ) = £ (lu c ) +ntu h + m (uj d ) . 
Using these expressions, the quasilinear potential (J73j) is written as 



(75) 



(76) 



(77) 



m£k 



Y 7T 5{u k - m- (u) m )) 



5K m£k cos 9 



(7S 



Note that wave-particle resonances involve harmonics of the bounce-averaged gyrofre- 
quency (u c ) and the bounce-averaged drift-frequency (u d ). By combining these results 
into Eq. (}4*S|) . we finally obtain the expression for the relativistic quasilinear diffusion ten- 
sor 



D 



QL 



E 

m,fc 



, ) 

i m£ muj k m 2 ) 



( £ 2 £uj k £m 
uj k £ (Jf. u k m 



7T 5(u k - m- (w m )) 



SK m£k cos 9 



(79) 
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The structure of the relativistic quasilinear diffusion tensor (|79|) clearly shows that, for 
arbitrary values of gyro-harmonic and drift-harmonic numbers (£, m) and non- vanishing 
wave frequency u>k, off-diagonal quasilinear transport cannot be neglected (see Appendix 
A for a brief discussion of off-diagonal quasilinear transport coefficients). Note also that 
wave-particle resonances, in fact, involve bounce-averaged cyclotron and drift frequencies, 
not their local expressions as might be expected. 

V. SUMMARY 

The present work has presented a complete derivation of the relativistic bounce-averaged 
quasilinear diffusion equation associated with arbitrary-frequency electromagnetic fluctua- 
tions in axisymmetric geometry. The main results are the relativistic quasilinear equation, 
Eq. (JTTj) and the corresponding diffusion tensor, Eq. (f79~|) . where SK mik is given by Eq. (ffiljj) 
or Eq. (jnZJ). In Eq. (f7U|) the s-dependent angle 9(s) is given by Eq. (J72J), (...) denotes 
bounce-averaging, and the Fourier-Bessel components (denoted by a tilde) are defined in 
Eq. (JoTj) . The quasilinear diffusion tensor is also presented in energy, pitch-angle coordi- 
nates in Appendix A. 

Future work will include a comparative study with previous analytical models of pitch- 
angle and energy diffusion 0|B], and with calculations of radial diffusion coefficients [2112110], 
plus generalization to allow non-axisymmetric unperturbed electric and magnetic fields. In 
related work, a numerical implementation of the present relativistic quasilinear diffusion 
equation in a multi-dimensional simulation code is being developed (in collaboration with 
Dr. Jay Albert), including an investigation of the importance of off-diagonal quasilinear 
diffusion coefficients in non-axisymmetric magnetic geometries. 
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Appendix A. Equatorial Pitch- Angle Formulation 

Since quasilinear transport is often discussed in terms of the equatorial pitch-angle 
instead of the gyro-action J g 2j, we introduce the cosine of the equatorial pitch angle 
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(denoted £) denned as 



£{J g ,e, Jd) 



JgUJ g0 



P 



e(7 + l)/2 



(a.i; 



where the parallel momentum pn and the rest-mass gyro-frequency u g o = qB Q (ip)/Mc are 
evaluated on the equatorial plane. From this definition, we obtain the differential relation 



where the drift frequency 



cO 



c^co dJg + uJdo dJd 



l-Hde 



J, 



du, 



cO 



defines an equatorial drift frequency [see Eq. (|B.2j) ] and 



PP||o 



go 



(A.2) 



(A.3) 



(A.4) 



is the Jacobian associated with the substitution J g — > £ (and thus has units of ac- 
tion). Using the differential relation (|A.2|) . the unperturbed guiding-center evolution of the 
equatorial pitch angle is £o = 0. 



The perturbed pitch-angle Hamilton's equation for 5£ is expressed as 
1 

d( dip 



5i 



d d , 2 x ( d d n 

^co — + u d0 — + (1 - 4 J 



The differential operator JF becomes 

'dF uj d0 dF \ d 1 dF d 



dt dt 



5H. 



T 



J^co d£ J dip J € <9£ d( 



de 



dF \ d 



J^ c0 <9£ J dt 



(A.5) 



Relativistic quasilinear diffusion equation (|4y|) can also be written in terms of the invariant 
coordinates I = (£, e, Jd) 



dFo 
dr 



1 8 (j<*K. bf » 



di J 



(A.6) 



where is defined in Eq. (jA.4|) and the components of the new quasilinear diffusion tensor 
are defined by the relation 



dl l nab df 



Qja QL QJb- 

The new quasilinear diffusion tensor is, therefore, expressed as 



D 



QL 



m,£,k 



( A, 



m£k 



UJk A, 



mlk 



mlk 



V m\ mek muj k 
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7™ Xmtk 
2 

m 



ralki 



(A.7) 



(A.8) 



where the pitch-angle coefficients A m g k and \ m e± ar e defined as 



^mtk 



I 2 , i 2 



cO 



£w c0 - (1 - £ 2 ) uj k + muj d0 



(A.9) 
(A.10) 



Note that, while the diagonal coefficient (|A.9J) cannot vanish (since it is strictly positive), 
the off-diagonal coefficient (jA.10|) may be small when evaluated at the wave-particle reso- 
nance u k = £(u) c ) +nu^+m (o?d). Further discussion of the importance of these off-diagonal 
quasilinear transport coefficients is, however, outside the scope of the present work and will 
be a subject for future work when non-axisymmetric magnetic geometries are also consid- 
ered. 



Appendix B. Bounce- Averaged Drift Frequency 

In this Appendix, we derive an explicit expression for the bounce-averaged drift fre- 
quency based on the definition (fTTj) : 

r ha' Vhl 9J d <l\ V \\\ 9s q OS J 

Using the identity 

c hl 9a = d_ (c \ _ ca d\p^ = 3_ (c \ + cJ^ du^ 
q ds ds \q J q ds ds \q J q\v\\\ ds 

which follows from the definition (jHJ) for p», we find that the second and third terms in 
Eq. (|B.1|) cancel each other and, hence, only the first term in Eq. (|B.1|) remains. The 
bounce-averaged drift frequency is, therefore, given as 



Based on this expression, we defined in Eq. flA.3|) the equatorial drift frequency cuao- 
An alternative expression for (u&) is obtained by introducing the bounce action 



and the bounce frequency 

r h V 9e 



= — = ^- , (B.4) 
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so that the bounce-averaged drift frequency (|B.2|) can also be expressed as 
where Eq. (JHJ) was used for p\\. 
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